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A complication in proving factorization theorems in Feynman gauge is that individual graphs 
give a super-leading power of the hard scale when all the gluons inducing the hard scattering are 
longitudinally polarized. With the aid of an example in gluon-mediated deep-inelastic scattering, 
we show that, although the super-leading terms cancel after a sum over graphs, there is a residual 
non-zero leading term from longitudinally polarized gluons. This is due to the non-zero transverse 
momenta of the gluons in the target. The non-cancellation, due to the non-Abelian property of 
the gauge group, is necessary to obtain the correct form of the gluon distribution function as a 
gauge-invariant matrix element. 
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I. INTRODUCTION 

Many studies of the factorization theorems of per- 
turbativc quantum chromodynamics (pQCD) were per- 
formed in axial or light-cone gauge. These gauges are 
appealing because there are no unphysical gluon polariza- 
tions or Faddeev-Popov ghosts, and the number density 
interpretation of parton distribution functions (PDFs) 
is clearer in light-cone gauge than in covariant gauges. 
However, in axial or light-cone gauge, the gauge-field 
propagator has unphysical singularities, which obstruct 
contour deformations needed in the proofs of factoriza- 
tion. Thus it is necessary to examine factorization in 
Feynman gauge where analytic properties associated with 
relativistic causality are manifest. A classic example of 
this is in the factorization proof for the Drell-Yan pro- 
cess 01 • More recently, the process-dependent Wil- 
son line directions in spin-dependent processes have been 
shown to be important for obtaining the correct relative 
sign in single-spin asymmetries [3, Q ; the derivation of 
the Wilson line directions is easiest in Feynman gauge. 

An essential complication in Feynman gauge is that 
exchanges of longitudinally polarized gluons between a 
collinear subgraph and the hard-scattering subgraph con- 
tribute at one higher power of the hard scale than in ax- 
ial gauge. One consequence is that there is no suppres- 
sion for adding extra gluon exchanges between the tar- 
get and hard-scattering subgraphs. A second and severe 
consequence is that when a hard scattering is induced 
by gluons, we find "super-leading" contributions, graph- 
by-graph. For example, super-leading terms for the deep 
inelastic structure function F2 are of order Q^, compared 
with the normal scaling behavior of Q° (up to logarithms 
of Q). Arguments using gauge invariance are needed to 
show that the super-leading terms cancel after a sum 
over graphs for the hard scattering, and that the leading 
terms combine to form gauge-invariant parton densities 
times conventional hard-scattering coefficients. Unfor- 
tunately the details are not clearly worked out in the 
literature. That there is a cancellation of super-leading 
terms was shown by Labastida and Sterman [5l|. They 



concentrated on issues in hadron-hadron scattering, but 
their argument applies more generally, for example to 
deep- inelastic lepton scattering (DIS). 

However, their argument also appears to suggest that 
there is a vanishing of the sum over graphs in which all 
exchanged gluons are longitudinally polarized. In fact, as 
we will show in this paper, the sum is actually nonzero; 
only the super-leading part cancels. This is a specific 
property of a non-Abelian gauge theory; the sum is in- 
deed zero in an Abelian theory. This is important be- 
cause the non-zero sum combines with terms with other 
polarizations for the gluons to give the correct gauge- 
invariant form Q\ for the gluon density. 

As far as we know, there is no detailed treatment in 
the literature of the role played by the gluon polarization 
in factorization, even at the level of low order Feynman 
graphs. In this paper, we give a detailed discussion at the 
level of two-gluon exchange (in the amplitude). Under- 
standing the exact nature of the contributions from the 
various kinds of polarization is important, not only to 
make sure that ordinary factorization is understood, but 
also to allow correct generalizations of factorization to be 
obtained (to include unintegrated PDFs, for instance). 

The key observation is that the target gluons are only 
approximately collinear — they have small but non- 
vanishing transverse components of momentum relative 
to their parent hadron. In contrast, the Labastida- 
Sterman argument treated the gluons as having zero 
transverse momentum. This was sufficient to show can- 
cellation of the super-leading terms. But it does not al- 
low a correct treatment of the remainder, which is power 
suppressed relative to the super-leading terms, but does 
contribute to the leading power. 

It should be emphasized that these issues are of direct 
practical importance. Perturbative calculations often re- 
quire the accurate identification of a gluon PDF. In a sub- 
tractive formalism, perturbative corrections to the PDF 
are needed to obtain correct higher order hard-scattering 
coefficients. In different formalisms it is important to 
ensure that the same gluon PDF is being used. In the 
pQCD dipole picture, for example, the qq dipole cross 
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section depends on the gluon PDF in the target (see Q 
and references therein). Ideally, it should be possible to 
show that this is the same gluon PDF that arises in stan- 
dard pQCD factorization theorems. 

The basic steps for deriving factorization for inclusive 
7*p — > X scattering are as follows: 

1. Consider the most general type of region that con- 
tributes at or above the standard leading power in 
Q/K. Here, A is a typical hadronic mass. In this 
paper, we will only consider the contribution illus- 
trated by the cut diagram in Fig. [Ija) . The hard 
scattering is induced by gluons from the lower bub- 
ble, £, all of whose lines are collinear to the tar- 
get. The final state of the hard scattering consists 
of a quark and antiquark, which emerge from the 
hard-scattering bubbles, Ti-L/R at wide angles. Ar- 
bitrarily many target collinear gluons (represented 
by the gluons and the "• • • ") may attach the target 
bubble to the hard bubbles. In general, the out- 
going quark and anti-quark evolve into bubbles of 
final-state collinear lines. However, the integration 
contour in the sum over final-state momentum may 
be deformed into the complex plane to where the 
final-state lines may be treated as being off-shell by 
order |5|. 

2. Use Ward identities and appropriate approxima- 
tions to disentangle and factorize all the "extra" 
gluons. For the set of graphs we consider, the re- 
sult should be a convolution product of the on-shell 
amplitude for qq scattering with a non- 
perturbative factor containing target-collinear glu- 
ons. A graphical representation of the result of this 
last step is shown in Fig.[TJb). The double lines rep- 
resent Wilson lines that make the parton densities 
gauge invariant. 

3. The non-perturbative factor in the final factoriza- 
tion formula of step 2 indicates an appropriate op- 
erator definition for the gluon PDF. This same op- 
erator definition appears in other factorization for- 
mulas for other processes, thus allowing the gluon 
density to be fitted and used in multiple calcula- 
tions. 

We will follow these steps for the case of one and two 
gluons. In a complete derivation of factorization, we also 
need to perform double counting subtractions (not shown 
explicitly in Fig. [T]). However, the details of the sub- 
traction procedure are not important for the main issues 
discussed in this paper, which focuses only on the factor- 
ization and identification of the gluon PDF for a specific 
region of momentum space. We will therefore include 
only a brief discussion of subtractions, in an appendix. 

In Sect, ini we describe the kinematics of DIS and the 
notation and conventions to be used in the rest of the 
paper. In Sect. IIIII we treat the simple case of a single 
target-collinear gluon, to establish our basic technique. 



In Sect. llVl we describe the properties of the standard in- 
tegrated gluon distribution function. Then in Sect.lVlwe 
extend the argument of Sect.|III]to the case of two target- 
collinear gluons and illustrate the role the longitudinally 
polarized gluons play in the definition of the gluon PDF 
in Feynman gauge. We give comments and a summary 
in Sect. IVII In an appendix we review the subtraction 
procedure appropriate for a more general treatment. 



II. BASIC SETUP 

The basic kinematic variables are the standard ones 
for DIS. The light-front coordinates of the momenta of 
the incoming proton and photon are, respectively. 



p= I p+,i£,o* 



2P^ 



q = -xP^ 



2xP+' 



(1) 

Since we have kept the proton mass nonzero, the longitu- 
dinal momentum fraction, x, is not exactly equal to the 
usual Bjorken XBj , but is related to it by 



l + A/l + 4:g4j 



(2) 



We work in the Breit frame, where xP^ — Q/^/2. To 
characterize the forward direction, we define the exactly 
light-like vector 



nj = (0,1, Ot). 



(3) 



In this paper we restrict our attention to contributions 
from graphs of the form of Fig. [Ha) , and it will be most 
convenient to work at the level of the amplitude shown 
in Fig. [21 The upper part of the graph, denoted by U, 
describes the scattering of a virtual photon off gluons. 
The lower part of the graph, denoted by C, describes the 
emission of the gluons from the target proton. In general, 
the amplitude for scattering off N gluons in the target 
can be conveniently expressed as the contraction of U 
and C With the momentum labels shown in Fig. [51 we 
have 



(4) 



Here j runs from 1 to A^. We will use to label the elec- 
tromagnetic index, and fij to label the gluon indices. The 
quark and anti-quark momenta are li and I2 and emerge 
from the hard scattering at wide angles. The kj label the 
gluons connecting the two subgraphs. We will be exam- 
ining the contribution when each is in a neighborhood of 
the target collinear region: 



A2 



(5) 



We can treat both Z/^(^i, ^2; {%}) and C{P;{kj}) as be- 
ing obtained by a sum over all graphs with the requisite 
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FIG. 2: Basic amplitude with arbitrarily many target 
coUinear gluons. 




P' 

(a) (b) 
FIG. 3: Amplitudes for scattering off a single gluon. 



external lines. From the examination of examples, it can 
be seen that each graphical contribution to the amplitude 
is obtained A^! times. We compensate by defining 
£(P; {kj}) to include a factor l/N\. 

III. SIMPLEST CASE: A SINGLE TARGET 
GLUON 

We start with the simple case of the exchange of one 
gluon, as shown in the graphs of Fig. [31 The amplitude 
is of the form 

M(^i^)=W;,g^'^/;,,. (6) 

The subscript {Ig) denotes that one gluon is exchanged, 
and we have explicitly shown the numerator of the 
gluon propagator. To simplify notation, we will often 
omit the momentum arguments and the index v of the 
electromagnetic current. 

The gluon is coUinear to the target, and and I2 are 
at wide angle, so that the quark subgraph is the hard 
scattering. We will show that, after a sum over diagrams. 



the lower part of the graphs in Fig. [3] can be separated 
into a factor that can be identified as a contribution to 
the gluon distribution function. 

There are regions where the target gluon may go far off- 
shell, but in that case it should be regarded as belonging 
to the hard-scattering subgraph at one higher power of 
5j . For this paper, we restrict ourselves to the case where 
the hard-scattering coefficient is lowest order (LO) in g^. 

A. Power counting 

To understand the power counting in Q, for large Q/A, 
we first observe that the target is highly boosted in the 
plus direction, so that 

C+^Q/K, C-^K/Q, (7) 

up to an overall power of A. Since li and I2 are at wide 
angles, all components of hi are of comparable size, and 
given by the dimension of the quark part of the graph: 

U- r^Ut^Q^^ (8) 
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To obtain the structure function, we square the ampli- 
tude and integrate over the final states. The integral over 
the gluon momentum k is invariant under boosts from the 
target rest frame, and therefore does not contribute any 
power of Q. The Lorentz-invariant phase space integral 
for the qq pair. 



(27r)64|li||l2 



.{2Tr)^S'^^'>ih+l2-q + k)., 



(9) 



is dimensionless. Thus it also contributes a zero power 
of Q when the quark-antiquark pair is at a wide angle, 
where all components of li and I2 are of order Q. Thus 
the overall power of Q is that for the squared amplitude. 

From Eqs. (O, we see that the largest power 

of Q is from the term U~C~^ ^ Q/A, for individual 
graphs. Squaring this gives the previously mentioned 
super-leading power of in a structure function. In 
the next subsection, we will use a generalizable method 
to demonstrate cancellation of the super-leading terms in 
the amplitude. To give a correct treatment of the lead- 
ing terms (<5°), we need to work at next-to- leading power 
accuracy in Q relative to the super-leading power. We 
also need to show that the leading terms correspond to 
a correct definition of the gluon density. 



B. K and G gluons 

A convenient technique was introduced by Grammcr 
and Yennie in the context of QED Q. In this approach, 
the propagator numerator in Eq. Q is split into two 
parts: 



5" 



where 



k ■ nj 



k^' 

k ■ nj 



(10) 

(11) 
(12) 



In these definitions, we use the light-like direction nj 
which characterizes the outgoing jet direction. The ratio- 
nale for this choice is that the dominant term inU ■ C is 
the one obtained from the g'^ component of the gluon- 
propagator numerator. reproduces exactly this 

component, the remaining components of K give smaller 
contributions, and K has a factor fc^ , which allows gauge 
invariance to be applied. 



J 



We use Eq. pO|) to separate Eq. ([6]) into two parts 
which we call the "if-term" and the "G-term." We will 
also use the terminology of the exchange of "if-gluons" 
and "G-gluons". 

From Eq. ([5|), we determine the Q dependence that 
results from the individual components in the K-term: 



U-K+-L+ 
All others ' 



Q 
A' 

— , or smaller. 



(13) 
(14) 

(15) 



The term gives the super-leading contribution, as- 

sociated with the power of Q in £+. Changing to 
removes one power of Q, by changing fc+ to /c*. The re- 
sulting power, (5°, we identify as the standard leading 
behavior. 

As we will see explicitly, gauge invariance ensures that 
k dotted into the sum of all graphs in the upper bubble 
vanishes, so that 



graphs 



(16) 



Hence the if-term can be dropped, so that Eq. ([6]) gives 
M(i<,) =Z^^-G^'^/:^. (17) 

Only the terms with transverse components on fj! are un- 
suppressed in Eq. ([TT]). This is easily verified by defining 



(18) 



and recalling Eqs. ([7]). Checking each combination 
of indices we find. 



U-L+ = 0. 



(19) 
(20) 
(21) 



Therefore, the transverse term dominates, and the other 
terms are power suppressed or zero. 



C. Gauge-invariance calculation 

We now ^ive an explicit demonstration that 
Egraphs^M'^^ ^ - 0. The graphs of Fig. [3] give 



7 +7 



/i 



(22) 
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Then for the if-term we have 



k ■ nj 



7+7 



(23) 



to which we apply the following identities for ^: 

$ = -{t2-$-m) + {t2-m) in term 1, (24) 

= {l/i - ti- ra) + {t^+ m) in term 2. (25) 

Using the Dirac equation for u{l2) and v(li), we find that two terms vanish and the remaining terms exactly cancel, 



so 



= \-u{h)Yv{h)+u{h)Y^{h)] = 0. 



k ■ nj 



Therefore, only the G-term survives, so that 



(26) 



(27) 



We now have the power-counting of Eqs. p^ ~ (PT|) . It is important to recognize that we wrote the amplitude in 
the form of Eq. (j27p without making any approximations. It is very tempting to replace k by an exactly coUinear 
value at the start of the argument, both in the upper, hard-scattering subgraph, and in Eqs. (fTTj) and p2)) . The 
approximated momentum is on-shell and has zero transverse momentum. This replacement is particularly natural to 
make in the K term, and was made by Labastida and Sterman Q. However, then we would have made an error in 
each graph contributing to the if -term that is of order kt/Q relative to the super-leading power Q. That is, the error 
would be of the same order as the leading terms in the final result. Using this kinematic approximation in the K term 
also entails using it in the G term, and we would not have obtained the contribution to ^ involving the transverse 
components oi k. If we make the collinear approximation in the upper subgraph but not in K , then the cancellation 
of the K-ieims in Eq. ((23)) is no longer exact. (The error would again be of order kt/Q relative to the super- leading 
power.) These observations will be particularly important for obtaining the correct two-gluon contribution in Sect.fVl 



D. Leading-power Factorization 

The next step in obtaining the factorization formula for scattering off a single target gluon is to drop power- 
suppressed terms. This means we can drop the contribution from Eq. ()20|) and keep only the sum over transverse 
components in Eq. (j27p . Furthermore, since we are restricting to the region of k given by Eq. we can now 
substitute the approximate parton momentum, 

fc=(fc+,0,Ot), (28) 

into the upper bubble only. Thus the hard scattering is initiated, as is usual, by a massless on-shell parton of zero 
transverse momentum. The amplitude is then 

Mi^ig)=Y^W{hM-.k)C,{P-k)+o(^^ , (29) 

where the index j runs only over the two transverse components. In an unpolarized cross section or structure function, 
we square the amplitude and sum/integrate over final states. However, we have assumed above that the outgoing 
quarks are at wide angles, so that the collinear approximation (|28p for the gluon can be safely applied in the upper 
part of the graphs. 

But in the integral over all final states, there is also a leading contribution from the region where one of the 
intermediate quark lines is collinear to the target, i.e., (fc — li) ~ (fc+ — /j^, 0, Of) or (fc — I2) ~ {k'^ — I2 , 0, Oj). This 
contribution to the cross section is already taken into account at the parton-model level, and therefore an appropriate 
subtraction should be made to avoid double counting, according to the principles summarized in the Appendix. Thus 
the contribution of the graphs to the structure tensor takes the form 

W""""' r^J2 J dliW' ^ {hM^kW'- ^' Cj{P]k)C]\P]k)~ snhiY&ciiont^^^ . (30) 

3,3' 
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Here, dH denotes the complete integration measure including sums over final states. The exact form of the subtraction 
terms is found by making the parton-model approximation on the struck quark, appropriate for the region of phase 
space where one of the quarks is target collinear. (The detailed steps for making this approximation are reviewed 
in the introductory sections of Ref. [9|.) We denote the parton-model approximation on the upper part of the 
graph by T-py^U^ {li,l2',k)W '^{li,l2',k). Furthermore, in the sum over final states for the unpolarized cross section, 
W{li,l2] k)U^ ^(/i, Z2; k) is diagonal in j and /, so we may write Eq. (pO| as. 



Y,Cr{P;k)c\,{P-k) 



O 



A 



(31) 



where W = W{li^ h'/k). We note that the manipulations needed to obtain Eq. |3T|) can be made transparent by first 
writing the contractions of the upper and lower bubbles in terms of sums over transverse polarizations: 



J2i^^ihj2; k)CjiP; k) = Y,T.^'{h,h; k) {et,,Y (et.Y MP; fc), 
j hi s 

where we have introduced transverse polarization vectors, 



(et,,y = (0,1,0,0), (£*,,)' = (0,0,1,0). 



(32) 



(33) 



Then by the use of diagonality oi W {li,l2]k)W ''{li,l2',k) after the integral over the final state, the sums over 
polarizations produce Eq. ([3T|) . 

The first factor in the integrand of Eq. (|3ip is just the unpolarized squared amplitude for j*g — > qq scattering 
minus a subtraction term. For the rest of this paper, we will ignore the details of implementing the double counting 
subtraction since it only modifies the hard scattering but does not affect Ward identity arguments such as those of 
Sect- HirCI 

Since W {li,l2', k) depends only on the plus component of k after the approximation of Eq. (j3ip . the integration 
over k^ and depends only on the target bubble. Thus, we will ultimately identify the lowest order contribution to 
the gluon PDF with the factor, 



X 



dk-d'^kt 
~(2^ 



Y,^r{P;k)c],{P;k), 



(34) 



where represents a general sum/integral over the final states of the target bubble, and the factor k^ gives the 
standard normalization for the gluon density. The factorization formula that follows from Eq. (|3ip can thus be 
expressed graphically as 



W 



dk+ 
2k+ 




Cf (P; fc)£], (P; k) - subtraction terms + 0{A/Q). (35) 



The upper part of the graph now has the interpretation of an on-shell scattering amplitude for a virtual photon to 
scatter off a transversely polarized gluon with exactly collinear momentum. In a general treatment of factorization, 
the steps of this section need to be generalized to more than one gluon, and factors like Eq. ([M]) need to be shown 
to correspond to an operator definition of the gluon distribution function. Therefore, before we go further we must 
examine the general properties of the gluon distribution function. 



IV. THE GLUON DISTRIBUTION FUNCTION 



A. Operator Definition 



For Eq. (IBS)) to be a consistent factorization formula for the single-gluon case, Eq. ([M]) should correspond to 
the relevant contribution to the operator definition of the gluon PDF. The standard definition [6| of the integrated 
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momentum space gluon distribution function for the proton is 

2 



f,/,{^,P) ^Y.j 2^'"^'^"'' Q+nQ,w-,Q,)VAQ+'{Q) IP). (36) 

Here — k ■ nj/P ■ nj = / is the longitudinal momentum fraction of the proton carried by the target gluon. To 
make the definition exactly gauge invariant, we have inserted a Wilson line operator in the adjoint representation, 



Pa = P exp ~ig, 



^ dy-A+iO,y-,Ot)Tp^, (37) 



where the Tp (/3 — 1, 2, . . .8) are the generators of SU(3) in the adjoint representation and P is the path-ordering 
operation. The T/^ are related to the structure constants, //j^k, by (TJg)^^ — —ifp-yn- The color indices of Va are 
contracted with those of the gauge field-strength tensor, whose definition we recall: 

g^^'^iz) = d^A^^{z) - d^A^^{z) ~ g,!^p^A%{z)A^^{z). (38) 

Then examining how to apply Eq. (|36p in perturbative calculations yields Feynman rules for the gluon PDF. In 
perturbative calculations, it is convenient to substitute the following identity for the Wilson line operator, 



Pex^p(^-ig, dy-A+{0,y-,Ot)Tp^ 



Pexp ^ 









t 



Pexp -i.g, / dy'-A+iy'-)TA , (39) 







and to insert a complete sum over final states. Then Eq. ([55]) becomes 



/./p(^'^) = EEE / ^^^~'^''''"~iP\Sl'iO,^-,Ot) Pexp / dy- A+ {y- + )Tp 

V a 7 = 1 ^ I \ Jo 



\X)x 



X (X|Pexp {-igs £ dy'-A+{y'-)Ti^ ^+^ (0) |P). (40) 

The symbol 'Ylix i^icludes all sums and integrals over final states. Expanding the Wilson line in small coupling gives 
the relation, 

Pexp(-tgs dy-A+{y-)TA^l + Pj2^-^f^Y[ dyi A+{y-)Tp. (41) 



The Feynman rules for the gluon PDF are found by using Eqs. (|4ip . ([55)1 inside Eq. PO)) and directly applying the 
rules of ordinary perturbation theory. 

B. Lowest Order 

We now expand the operators in Eq. (|40p in powers of gs- At zeroth order, only the derivative terms in Eq. (|38p 
contribute to Eq. 

2 _ 

/9/p(C.n =EEE / ^^e-«^"-"(P|(5+A^Jz.-)-a^Ai(^-))|X)(X|(a+A^jO)-a^Ai(0))|P). (42) 

X a j=l J ^ 

Directly applying the rules of perturbation theory yields 
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where the factors and C"^ denote the parts of the incoming target bubble to the left and right of the final-state 
cut, with their color indices. Equation is exactly what we asserted to be the gluon-density factor, Eq. ([M)) . for 
DIS with a single target gluon. To express Eq. ((43|l diagrammatically, we define the following Feynman rules: There 
is a lower bubble representing the incoming target. 




(44) 



For the g^^ — k^n^/k ■ nj factors in Eq. (|43)) we use the notation 



k ■ rij 



(45) 



Since Eq. (|45p modifies the vertex to the gluon from the target bubble in Eq. ((43|l . we refer to it as the "special" 
vertex. Diagrammatically, we notate the contribution to the gluon PDF from a single gluon emission on each side of 
the cut as 




(46) 



Implicit in this notation are an overall factor of fP^, integrals of k and fc^, and a sum over the gluon color index a. 
Then Eq. ([35)l can be written as 






subtraction terms + O ( A/Q) . (47) 



C. Beyond Lowest Order 



General Feynman rules for the gluon PDF, Eq. (|36p. are found 6] by keeping other terms in the expansion of the 
Wilson line operator, and the third term in Eq. ([38|) . To start out, let us ignore the third term in Eq. ([38| . Looking 
only at the operators on one side of the final-state cut and using Eq. (|4T|) , we get factors in the integrand of the form 

- z ffci • - k{n't'\ '9.nTUi^2n, ^ ^9snf f ^ ^ ^ ^ ^ igsTi'','' f,^_,^^,^ (48) 

^ ' ^ YI!=2k^■n,-^Q Y^^,h■n,-^Q k^ ■ nj - zO 

At this order of perturbation theory, N gluons are emitted from the target bubble, with color indices ai, • • • , Un- 

Now we treat the contributions which use the third term in the field-strength tensor Eq. p8|) . To stay at the same 
order of perturbation theory as in (|48p . we must drop down one order in the expansion of the Wilson line operator. 
Thus, we obtain contributions of the form 
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FIG. 5: A general contribution to the gluon distribution function. 



Rather than treating the derivative terms and the third term in Eq. separately, we can combine (pS)) and (H^)) at 
each order of perturbation theory. The complete contribution is then 



The first factor in ([50| tells us what the special gluon vertex analogous to Eq. (j45|) is for N gluons. It differs from 
the corresponding factor in ((48|) by having the gluon momentum ki replaced by the total momentum k in the first 
term only. We therefore obtain the Feynman rules for the gluon PDF shown in Fig. 21 where factors —i and i on the 
left and the right of the final-state cut have been dropped, and we have removed a factor of — k ■ nj, just as in 
Eq. (|45l) . A general contribution to the gluon PDF is illustrated in Fig [3 To simplify expressions, in later sections 
we will drop the explicit appearance of — iO in the eikonal denominators. 



V. SCATTERING OFF TWO COLLINEAR GLUONS 



We now extend the treatment of Sect. Illll to the case of two target-coUinear gluons. The graphs under consideration 
are shown in Fig. [5) Our aim is to show that when the graphs are summed. Ward identities similar to what were 
used in Sect. Illll lead to a factorized structure. For a consistent factorization formula, the resulting gluon PDF should 
correspond to what is obtained using the Feynman rules of Fig. [5) 



A. Power counting and K, G decomposition 

As in the previous sections, we write the amplitude obtained from the graphs in Fig. [S] as the contraction of an 
upper part and a lower part, connected by gluon-propagator numerators: 

M^,gj^U,>^,y^'^^g'^'-^^-£^,^,. (51) 

All incoming gluons are target coUinear, with the sizes of the momentum components given in Eq. ([5]). The gluons 
have leading contributions from regions where they are far off-shell, but these contribute only at high order in g^. 
Furthermore, since there are only gluon-fermion couplings in the upper part of the graph, ghosts do not enter at this 
level of the calculation. A general, all-orders proof will need to take into account ghost attachments as well. 

The upper part has one more external gluon than before, so its dimension is reduced by a single power of Q; all the 
components of l^^'^^i'^ are therefore to be treated as order . The lower bubble now includes two Lorentz indices. 
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(d) (e) (f) 

FIG. 6: DIS with two target-coUinear gluons. 

and in the boosted frame has, up to overall powers of order A, the power dependence, 

A2' 
Q 
A' 

£" ~ £+- ^ C-+ ^ Q°, 
A 



A! 





C++ 




^C+' 


-+- ^ 










c— 



(52) 

(53) 
(54) 
(55) 

(56) 



Combined with the Q^^ factor from U, this gives a super-leading contribution from the ++ components, leading 
contributions from the +t and t+ components, and power-suppressed contributions from the other components. 

That the term gives a superleading result is an example of a general result from power-counting arguments: A 
superleading contribution arises whenever all the gluons joining the coUinear and hard subgraphs are longitudinally 
polarized. We will show that there is a cancellation after a sum over all graphs for U leaving a result that is just 
leading power. That there is a leading contribution from the t+ and +t terms is expected: They correspond to a hard 
scattering on a transversely polarized gluon accompanied by a longitudinally polarized gluon. The remaining terms 
are of a non-leading power, and we can therefore ignore them. We will show that after a sum over graphs, the total of 
the H — \-, +t and t-\- terms is equivalent at leading power to a hard scattering off a single transversely polarized gluon 
multiplied by the appropriate two-gluon factor obtained from the Feynman rules for the gluon PDF. 

Just as with the singlc-giuon case, we split the propagator numerators into K and G terms: 



g 



9 



where Ki, K2, Gi and G2 arc defined just like K, 
momentum, fci or fc2- So we write Eq. (|5ip as 



(57) 
(58) 



G, except for the replacement of k by the appropriate gluon 



(23) 



= 11, 



G^i 



(59) 
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Exactly as in the previous section, the K1K2 term gives a super-leading contribution, graph- by-graph, while the K1G2 
and G1K2 terms are leading. The G1G2 term is power suppressed, and so we neglect it. Each of the non-suppressed 
terms has a factor of Ki or K2 , so we will next apply gauge invariance to get a result of the general form 

M(2g) - Up, (fci + k2) gp'P {■■■ )p + (A/Q) . (60) 

The first factor is just the basic amplitude already considered in the previous section, but with gluon momentum 
ki +k2: 



Up,{ki + k2) = 




The remaining factor (• • ■) ^ will determine what "special" vertex, analogous to Eq. (|45p . we should use for the gluon 
density. 



B. Gauge invariance on K-z factor 



We first apply the method of Sec. IIII CI to a single i^T-gluon in Eq. ([5^ . The explicit expression for the upper 
bubbles in the first three graphs in Fig. [6] is 



+ 7 



Ml 



^2 



V^2 



jl 2 ^2 



-/i- 



vih). 
(62) 



We organize these graphs by having the ki gluon attached in one place on the quark line, and then the graphs 
correspond to the different possible placements of the /c2 gluon. Similarly, for the other three graphs. 



7 



1 



^1 + $2 /i 



vM2 



7' 



1 



1^2 — /l — ™ 



1 



/2 $2 ' 



vih). 
(63) 

In Eq. ([55]) . the terms that are leading or super- leading have at least one ii'-gluon, to which we apply the method 
of Sec. IIII Cl We start with the case of any term involving K2. For graphs (a)-(c), we have a factor 



(a-c) 



(2 - ^) "'^^ V/2 - 1^1 - 1^2 ~ 

7X^/3+ 7a.; (7 7 

\l2-fl-m 



1 



7Mi 



/2 



^2 ^1 ^' 



7 



/i 



^2^atf3 



k2 ■ nj 



vih). 



Calling the terms in brackets term 1, term 2, and term 3, we now make use of the following identities: 

^2 = ^{(2 ^ ^2 ^ ^) + (^2 ^ "f^) ij^ term 1, 

= —{(2 ^ ti ^ $2 ^ + {1-2 ^ ^ "^) term 2, 

= i$2 ^ ti ^ m) + (/i + w) in term 3. 
Using the Dirac equation to eliminate two of the terms, we then have 

1 



(64) 

(65) 
(66) 
(67) 



■ 7m i 



(a-c) 
1 



■ Qsuih) 



-7Mi 



(2 ^2 ^ 

1 



7 tpto 



^2 ~ 1^1 ~ ^2 ~ ^ \/2 ~ ^1 ~ ^ \/2 ~ ^1 



1 



vih)^. 
k2 ■ nj 



(68) 
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The last two terms cancel because of the Abelian nature of the QED coupling. However, the remaining terms lead to 
a non- vanishing result due to the non- vanishing commutation relations for and t^, 



(a-c) 



1 



Similarly for graphs (d)-(f) we have 



(d-f) 



Adding Eqs. ^ and we get 



AI2M2 



9sU{l2) 



^2 $1 ^2 



^2 + ^1 /l 



7+7 



igsfapK 



k2 ■ 



v{h) 



igsfa 



/3k 



(69) 



(70) 



(71) 



This is the general expression for the upper bubble when it is contracted with the i4r2-term from gluon 2. Notice that 
it vanishes when the theory is Abelian. This is an example of a general result that X-gluons give zero contribution 
in an Abelian gauge theory. Notice also that it is of the form of the upper factor ([221) for one gluon of momentum 
ki + k2, times an eikonal factor. 



C. Gauge invariance on Ki factor 



If we instead consider the contraction of the upper bubble with the Ki gluon from gluon 1, then exactly similar 
steps lead to, 



-gsuih 



7^2 K 



1 



^2 $2 ^ 



7 + 7 



1 



ii — fi — ■ 



7aI2 K 



v{h) 



igsfafSK 



ki ■ nj 



.(72) 

Note the overall minus sign that arises from the reversed roles of gluon 1 and gluon 2. To complete the analysis, we 
must consider the separate cases: where the other gluon is also a iC-gluon, or where one gluon is a X-gluon and the 
other is a G-gluon. 



D. K1K2 term 

The K1K2 term in Eq. (|59p follows immediately from Eq. (j7ip when we include the contraction with Ki: 

1 \ „ „ / 1 



gsuih 



12 



— — m 



7 +7 



^2 ^ - i\ — m 

^ '^gs faPnki 



v{h)x 



M "J 



(fci • nj){k2 ■ nj) 



(73) 



This has the form of the amplitude (|22p with a single gluon of momentum fci + k2 multiplied by a special vertex. 
Thus for the K1K2 term in Eq. ([55)) . we have: 



^;f,i^r'^"'^2^"^^Mf.2 = + ^2) 



igsfaPnkl.p 



'J "J 



(fci • nj)(fc2 • nj) 



C 



a/3 



(74) 



The two ii'-gluons have been factored out into eikonal couplings. A diagrammatic representation of Eq. (|73p is shown 
in Fig. [71 Note that this term would exactly vanish in an Abelian theory, where fa^K = 0. 

Graph-by-graph, this is still super-leading, since the factor of g^ ^ and everything to its right power-counts just 
like C for the case of one-gluon exchange. Furthermore, (|74p is not symmetric between the gluon momenta, fci and 
fc2, despite the symmetric occurrence of Ki and K2 factors on the left-hand side. Both problems are remedied by 
applying the K + G decomposition to the g'' ^ factor in Eq. ([7^ : 



y — -^(1+2) ^ ^(1+2)' 



(75) 
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FIG. 7: Graphical structure of the K1-K2 term. The arrows are used to represent Jf-gluons which couple to the quark via the 
special coupling in Eq. (|74|l . 



where 



K 



p p _ 

(1+2) 



(fcl 



"J 



nP P — nP P 

■-^(1+2) ~ y 



{ki+k2)P 

1 

(fci + fc2) • nj 



(76) 
(77) 



The K(i^2) term gives zero when multiplied with Up'{ki + (see the calculation in Sect. IIII C|) . This removes the 
super-leading part, just as in one-gluon exchange, leaving only the G(i+2) term. Hence, the K1K2 term is 



Z^p'(fci+A:2)G^i;2) 



iQsfaPuki^p - 



(fc2 • nj)(fc2 • nj) 
= Upi{ki + ^2^ 

Thus it has the structure of Eq. ([SO)) . 



i9sfaf3KnfnP'' 



(fci • rij)(fc2 • »^j)(fci + fe) • nj 



(fcf'(fc2-n,)-fc2^'(fci.nj))/:^f^,^. (78) 



E. G1K2 term 



We next apply Eq. ([7T|) to the G1K2 term to obtain 



QfJ'lPl j^P2P2 jf^aP 
^1 ^2 ^ MIM 



--U^{ki+k2) 



igsfaf3K 9 



,P2 



kPnP' 
ki -nj J k2 ■ nj 



/•a/3 



(79) 



again with the structure of Eq. (|60p . As already observed, there is no super-leading contribution from the G1K2 term, 
only a leading-power contribution. 



F. ^^102 term 



For the K1G2 term we similarly obtain 



k2-n,jj fci • nj 



a/3 



(80) 



Note the overall minus sign relative to Eq. (I79p that arises because of the reversed role of gluons 1 and 2. A graphical 
representation of the factorization in Eqs. (|79p and ([50)1 is shown in Fig. [51 



G. Factorization and the Gluon Distribution Function 



Adding Eqs. ([28|), ([79|), and ([80|) gives 



M(2g) =z^;(fci 



igsf, 



a(3K 



J 



^2 • fci ■ (fci -I- fc2) • nj \fci • nj fc2 • rij 



^2 



kP 



pa/3 



0(A/g). (81) 
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(a) 





(b) 

FIG. 8; Graphical structure of (a) the G1-K2 term and (b) the G2-K1 term. 



This may be written compactly if we define the factor in braces to be 



C 



- r 

(2g) ~ ''■9sJapK 



k2 ■ nj ki ■ nj (fci + k2) ■ nj \ ki ■ nj k2 ■ nj 



k^ 



K 



c 



a/3 



(82) 



Normally a vector like this would have its p = + component power counting as Q/A. In fact this component is zero. 
The reason is that it is constructed from a combination of Gi, 62, and G(i_|_2), each of which gives this property 
individually. 

The importance of the last remark is that it shows that the dominant contribution, i.e., the leading-power contri- 
bution, is from transverse values of p, so that lA^ which goes into the hard-scattering factor, can be treated as having 
an incoming transversely polarized gluon. Compare this with what is found in Eq. (j29[) . 



It is the factor i^'^^^g) which we would like to identify with a two-gluon factor in the amplitude for the gluon 
distribution function. The derivation so far, where we have systematically eliminated super-leading contributions, 
shows that we have a result with the standard leading power of Q, i.e., . The leading terms are for p = j, a 
transverse index: 



'^j(23)(-P;'=l'^2) = igsfafiK 



k 



k2-nj ki ■ rij ' (fci + fc2) ■ V^i ' "-J ^2 ■ rij 



2j 



kij 



(83) 



The steps for obtaining a factorization formula are now exactly analogous to the steps in Sect. IIIIDI The analogue 
of Eq. ((29)) for the case of two gluons is 



Af(2<,) =Y,W{hM: k)^2g){P-^ k2) 
3 = 1 



ki + k2 





-OiA/Q). (84) 



Here we have defined k = ki + k2. 



The factor ^''(^2g) ^ii ^2), given in Eq. is exactly 
what is obtained from the Feynman rules for the gluon 
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PDF listed in Fig. [4l This would not have happened if 
we had neglected the K1K2 contribution from Sect. IV Dl 
to the leading behavior in Eq. ((83|) . Equation (|78|) is 
therefore important for agreement between the Feynman 
rules for the giuon density and the leading contribution 
to DIS. 

To summarize this section, we find that leading-power 
contributions arise when one or both gluons is longi- 
tudinally polarized. With one longitudinally polarized 
gluon, we get expected eikonal factors in Eqs. (|79l80p . 
after a sum over graphs. These factors are analogous 
to the eikonal propagators that appear in the quark 
PDF, but they do not have the correct Feynman rules 
to correspond to the gluon PDF. A further contribution 
Eq. ([78]) arises when both gluons are longitudinally po- 
larized. Previously, this contribution was expected to 
vanish (see, e.g., Q). But it is in fact non- vanishing, and 
is needed for a correct correspondence with the definition 
of the gluon PDF. 



VI. SUMMARY AND CONCLUSION 

We have given a direct illustration using gluon-induced 
DIS that contributions from longitudinally polarized glu- 
ons do not cancel in sums over graphs in Feynman gauge. 
Indeed, they yield leading contributions that are needed 
to correctly identify the standard gluon PDF. Calcula- 
tions that go beyond lowest order in gs therefore require 
care in how the polarization of external gluon lines is 
treated, and involves keeping appropriate combinations 
of K- and G-terms in a Grammer-Yennie style treatment. 
Although we have obtained our result for the simple case 
of the integrated PDF, our arguments relied only on the 
application of Ward identities at the amplitude level. 
Therefore, similar results should hold for unintegrated 
PDFs (and related objects like fragmentation functions) 
once suitable definitions have been established. Since 
Eq. ([71)1 exactly vanishes in an Abelian theory, this re- 
sult is a specific example of how the non- Abelian nature 
of QCD can complicate factorization arguments. 

Higher-order calculations of hard-scattering coeffi- 
cients involve subtraction terms with multiple target- 
collinear gluons, so the results obtained here are of direct 
importance for obtaining correct hard-scattering coeffi- 
cients beyond lowest order consistent with factorization 
and a well-defined gluon PDF. 

Standard power-counting shows that it is also possi- 
ble to get leading contributions when the hard scatter- 
ing is induced by Faddeev-Popov ghosts, but only at 
higher order in the hard scattering than is treated in 
this paper. It is natural to expect a generalization of 
the known results [13, lUl for the operator product ex- 
pansion that applies to moments of DIS structure func- 
tions. On this basis we expect the general result to be 
that the sum over all graphs gives extra contributions to 
the factorization properties involving coefficients times 
what we will call alien PDFs. The operators defining 



/ 

FIG. 9: Handbag diagram for DIS. 




FIG. 10: Parton model approximation to handbag diagram. 

the alien PDFs are variations of certain parent oper- 
ators under Becchi-Rouet-Stora-Tyutin transformations 
and therefore vanish in physical matrix elements. Non- 
vanishing of the Green functions of the alien operators 
with off-shell quark and gluon states has caused calcula- 
tional problems [l^, [l^ . 
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APPENDIX A: SUBTRACTIONS 

We now review the subtractive approach — e.g., [Tsj 
and Sect. VI of 0] — we use for determining the hard- 
scattering factorization. It is generalized from the Bo- 
goliubov approach to renormalization. Our strategy is 
to recursively examine successively larger leading regions 
for the DIS cross sections. 

The procedure starts with a handbag-diagram struc- 
ture — Fig. [9] — to give the LO term, i.e., the parton- 
model formula. This term is obtained by an approxima- 
tion valid to leading power when the incoming quark is 
nearly on-shell and coUinear to the target. The approx- 
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FIG. 11: Generalized handbag diagram, with arbitrarily many 
extra gluon exchanges. 



imation is denoted by the hooked Une in Fig. \W[ and it 
involves replacing the momentum of the struck quark in 
the upper part of the graph by its parton-model approxi- 
mation, which is massless, on-shell, and of zero transverse 
momentum. 

If a general graph contributing to the cross section is 
denoted by F, then we call graphs where we make the 
parton-model approximation TloT- At this level we also 
need to consider graphs with arbitrarily many target- 
coUinear gluons attaching to the the hard vertex — Fig. 
[TT] A Ward identity allows the target-coUinear gluons to 
be disentangled from the LO hard-scattering coefficient. 
Since there are no super-leading contributions in these 
graphs, and since the relevant graphs are tree graphs, 
the Ward identities are unproblematic. 

The result is a convolution product of the LO hard- 
scattering coefficient with a sum of graphs identifiable as 
an expansion of the quark PDF. Schematically, the cross 
section is then written, 



A\" 



a = TLo5]r + 0(g» + o((-j a) = 



A 



(Al) 



where fq/p{x, Q^) is the quark PDF, Clo is the LO hard- 
scattering coefficient (from the electromagnetic vertex), 
and ® symbolizes the usual convolution product. 

Some of the errors in this approximation are power 
suppressed, indicated by the term [K/QY where a > 0, 
and we do not consider them further. Other errors are 
caused by regions with larger transverse momentum for 
the quarks and gluons, and by graphs not of the form of 
Figs. [H] and [TT] These will be covered by our treatment 
of higher-order scattering, and are suppressed only by a 
power of the strong coupling at scale Q, as indicated by 
the term O(g^cr). 

Hence, we have the LO approximation to the DIS cross 
section 



CTLO = Clo ® fq/p(x, Q^)- 



(A2) 



To find the next-to-leading (NLO) contribution, we ex- 
amine, among others, graphs of the topologies shown in 






(a) 



FIG. 12: Examples of graphs giving NLO contributions for 
DIS. These are the relevant ones for the gluon-induced term. 
But note that there are also other graphs that give the quark- 
induced NLO contributions. 




A- 





-4 











FIG. 



13: Generalized parton-model approximation to Fig. 
Each graph corresponds to a term in Eq. (|A4|I . 



Fig. [TH Specific examples are included in the graphs we 
treat in the main body of the paper. For the moment 
we ignore the possibility of extra gluon exchanges. Some 
of the graphs, e.g.. Fig. [T^ a). are among those already 
considered in the LO approximation, merely viewed with 
a different partition between the lower bubble and the 
parton lines explicitly drawn at the top. 

Therefore, before we can get the true NLO contribu- 
tion, we must first subtract the LO approximation: 



(A3) 



This is just the remainder from applying the approxima- 
tion already considered in Eq. (|Aip . Finally, to evaluate 
Eq. (|A3|) , we apply approximations that are good for the 
wide-angle 2^2 parton subprocesses, denoted by the 
symbol Tnlq: 

TnLO (E r - ^LO E r) = E ^NLO,, ® fj/piQ^), 

(A4) 

where we must now allow for a sum over parton flavors. 
For the graph in Fig. [TWa) , the result is shown graphi- 
cally in Fig. [131 

If graphs like Fig. [T^ a) were all that mattered, then 
this would complete the discussion of the subtraction for- 
malism. However, the gauge-invariant LO term needs a 
generalized handbag structure including extra gluon ex- 
changes, as in Fig. 1111 There are thus corresponding 
subtractions for graphs like Fig. I12f b). as shown in Fig. 
1141 Implementing the appropriate approximations and 
applying Ward identities allows us to identify how the 
second term in Fig. [T^ contributes at LO. After applying 
Ward identities, the second term in Fig. [Hjcan be written 
diagrammatically as in Fig. 1151 There we see the sepa- 
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ration of the graph into the LO hard-scattering matrix 
element and a single-gluon contribution to the PDF. 



FIG. 14: NLO term corresponding to Fig. [T2T b). with sub- 
traction for non-handbag LO contribution. The subtraction 
term uses the same graph as the first term, but it has been 
drawn differently to show that it gives a case of Fig. 1111 




FIG. 15: Particular graph for gluon-induced NLO term, with 
subtraction for non-handbag LO contribution. 



Thus we see that evaluating Eq. (|A4[) for the NLO cor- 
rection requires that we know exactly the single-gluon 
correction to the quark PDF in Eq. (|A2|) . since the defi- 
nition of the approximation Tlo ^ is what gives us the 
lowest order, parton- model factorization in Eq. (|Aip . In 
other words, using the subtraction approach to calculate 
higher-order corrections to the hard-scattering coefficient 
requires that we know exactly what we are subtracting. 



[1] J. C. Collins, D. E. Soper and G. Sterman, Nucl. Phys. B 
261, 104 (1985); J. C. Collins, D. E. Soper and G. Ster- 
man, Nucl. Phys. B 308, 833 (1988). 

[2] G. T. Bodwin, Phys. Rev. D 31, 2616 (1985) [Erratum- 
ibid. D 34, 3932 (1986)]. 

[3] J. C. Co llins, Phys. Lett. B 536, 43 (2002) 
[arXiv:hep-ph/0 204004^ . 

[4] S. J. Brodsky, D. S. Hwang and I. Schmi dt, Phys. Lett. 
B 530, 99 (2002) 'arXiv:hep-ph/0201296\ 

[5] J. M. F. Labastida and G. Sterman, Nucl. Phys. B 254, 
425 (1985). 

[6] J. C. CoUins and D. E. Soper, Nucl. Phys. B 194, 445 
(1982). 

[7] C. Ewerz, A. von Manteuffel and O. Nachtmann, Phys. 

Rev. D 77, 074022 (2008) [arXiv:0 708.3455 [hep-ph]]. 
[8] G.J. Grammer and D.R. Yennie, Phys. Rev. D 8, 4332 



(1973). 

[9] J. C. Collins, T. C. Rogers and A. M. Stasto, Phys. Rev. 
D 77, 085009 (2008) arXiv:0708.2833 [hep-ph]]. 
[10] S. D. Joglekar and B. W. Lee, Annals Phys. 97, 160 
(1976). 

[11] M. Henneaux, Phys. Lett. B 313, 35 (1993) [Erratum- 

ibid. B 316, 633 (1993)] arXiv:hep-th/930610L . 
[12] J. C. Collins and R. J. Scalise, Phys. Rev. D 50, 4117 

(1994) arXiv:hep-ph/9403231 . 
[13] B. W. Harris and J. Smith, Phys. Rev. D 51 4550 (1995) 

arXiv:hep-ph/9409405 . 
[14] D. Binosi and L. Theussl, Comput. Phys. Commun. 161, 

76 (2004) arXiv:hep-ph/0309015 . 
[15] J.C. Collins, Phys. Rev. D 65, 094016 (2002) 

[arXiv:hep-ph/0110l"T3] . 



